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OBSERVATIONS  OF  ELECTRON  TEMPERATURE  GRADIENTS 

IN  MID-LATITUDE  E -LAYERS 

s 


I.  INTRODUCTION 


The  underlying,  mechanism  for  the  formation  of  mid- 
latitude sporadic-E  layers  by  the  convergence  of  metallic 
ions  to  wind  shear  nodes  (Axford,  1967;  Macleod,  1966 • 
Whitehead,  1967)  is  found  to  account  for  much  of  the 
observed  data  (Matshushita  and  Smith,  1972  and  1975),  How- 
ever, the  temperature  predictions  of  individual  particle 
populations  in  and  near  the  layers  (Hines,  1965-  Gleeson 
and  Axford,  1967;  Hooke.  1969a, b)  have  not  been  fully 
tested,  largely  as  a result  of  inadequacies  in  diagnostic 
devices  or  complementarity  of  payload  instrumentation. 

The  work  of  Gleeson  and  Axford  (1967),  and  Hooke 
(1969)  describe  the  processes  that  determine  charged- 
particle  temperatures  at  E-region  altitudes,  and  they 
employed  their  results  to  estimate  the  influences  of 
electron  density  and  ion  composition  on  height  variations 
of  plasma  temperatures  in  mid-latitude  sporadic-E  layers. 
They  predict  minima  in  both  T and  (T  - T ) within  the 
layer,  with  bottomside  values  substantially  lower  than  the 
topside  (Hooke,  1969a).  The  two  cases  described  in  detail 
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by  Hooke  show  that  topside  electron  temperature  can  exceed 
that  at  the  layer  peak  by  at  least  150”k  whereas  the 
electrons  on  the  bottomside  can  be  isothermal  (within  lO'^K)  . 
In  the  work  of  Hooke,  allowances  were  made  for  variations 
in  the  neutral  Ras  temperature  that  can  accompany  tidal 
and  Rravity-wave  motions.  He  clearly  shows  that  variations 
in  can  play  a major  role  in  the  final  height  variation 
of  electron  and  ion  temperature  through  layers. 

In  contrast  to  the  works  of  Gleeson,  Axford,  and 

Hooke,  which  allow  T ^ T , is  the  position  that  all  parti- 

cles  below  130  km  are  in  thermodynamic  equilibrium,  where 

variation  in  the  temperature  of  the  neutral  species  is 

accompanied  by  identical  variations  in  the  temperature 

of  the  charged  particle  population.  To  couple  this  concept 

with  a proposed  mechanism  of  reversible  heating  by  gravity 

waves  and  tides  (Hines.  1965'  Hooke  1969  a.b)  requires  that 

changes  in  T cross  an  E layer  be  identified  with  the  T 
e s n 

variation  produced  by  the  gravity  waves. 

In  this  work  we  present  pu Ised-plasma-probe  measure- 
ment.s  of  electron  temperature  which  show  that  substantial 
temperature  gradients  do  exist  within  an  E^  -layer.  This 
data  was  collected  aboard  an  Aerobee  150  rocket  launched 
from  White  Sands  Missile  Range  on  29  June  1974  at  2031  MST, 
The  payload  reached  an  altitude  of  229.5  km  and  up-  and 
downleg  portions  of  the  flight  passed  through  blanketing 
sporadic-E  layers  centered  at  106.7  and  106,1  km. 
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respectively.  The  layers  are  shown  in  Figure  1 where  the 

flight  profiles  of  relative  electron  density  are  measured 

by  peak  electron  currents  collected  by  the  probe.  The  up- 

and  downleg  E thicknesses  at  the  0.1  point  were  1.16 

s e 

and  .73  km  with  corresponding  peak  electron  densities  of 
3(10^)  and  7(10^)  cm  respectively. 

Within  the  layers,  the  probe  simultaneously  measur- 
ed changes  in  ambient  plasma  density  and  electron  tempera- 
ture with  a spatial  resolution  of  14  and  300  meters, 
respectively.  This  capability  made  temperature  measure- 
ments possible  within  tlie  steep  density  gradients  char- 
acteristic of  sporadic-E.  This  paper  concentrates  on 
these  temperature  results,  showing  that  the  peak  and 
bottomside  values  of  T^  were  approximately  equal  at  340°K 
but  with  the  topsides  showing  T^  enhancements  with  gradi- 
ents as  liigh  as  34o”k  km.  Before  discussing  the  details 
of  these  measurements  and  their  implications  with  regard 
to  theoretical  concepts,  we  describe  the  experimental 
technique  and  analysis  procedure. 

II.  E.XPER  I MENTAL  TECHNIQUE 
The  instrument.  A pu Ised-plasma-probe  (referred  to 

3 

as  P ) was  utilized  for  the  simultaneous  measurements  of 
electron  density,  temperature,  and  density  variations.  The 

3 

P is  a Langmuir-type  probe  with  a specialized  electronic 
procedure  for  generating  current-voltage  characteristics 
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in  a way  which  results  in  improved  reliability  and  expanded 
versatility  in  laboratory,  ionospheric  and  reentry  plasma 
studies  (Holmes  and  Szuszczewicz , 1975;  Szuszczewicz  and 
Holmes,  1975;  1976a, b).  Compared  with  conventional 

3 

Langmuir  probe  procedures,  the  P technique  reduces  sub- 
stantially the  errors  that  can  affect  determinations  of 
electron  densities  and  energy  distribution  functions  when 
probe  surface  conditions  change  within  the  measurement 
period  (Hirao  and  Oyama , 1972;  Sturges,  1973;  Szuszczewicz 

3 

and  Holmes,  1975;  Wehner  and  Medicus,  1952),  P diag- 
nostic capabilities  extend  into  the  domain  of  temporal  and 
spatial  plasma  variations  where  the  technique  additionally 
provides  a simultaneous  measurement  of  density  fluctuation 
I)ower  spectrum  and  demonstrates  itself  as  an  important 

flevice  for  in  situ  investigations  of  irregular  ionospheric 
s t rue  t ures . 

The  probe  was  a straight  16.4  cm  length  x 0.06  cm 
diameter  tungsten  wire  extended  from  the  payload  tip  on 
a 57  cm  long  boom  electrically  tied  to  the  rocket  potential. 
Tfie  probe  was  separated  from  the  boom  by  a cylindrical 
guard  electrode  of  1,6  cm  length,  and  both  the  probe  and 
ttie  guard  electrode  were  swept  with  identical  electric 
potentials  of  the  form  shown  in  Fig.  2 and  specified  by  the 
parameters 
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380  msec 


(T  T,J 

oil,  M 


(100ms,  1900  /is) 


(V_,  V^)  = (-2.38V,  + 3.02v) 


V,J  - 0 V 

As  shown  ill  the  figure,  two  alternating  modes  of 
probe  operation  were  employed.  Figure  2A  depicts  a linear 
sawtooth  voltage  which  represents  the  conventional  approach 
to  Langmuir  probe  operation  while  Fig.  2B  shows  the  pulse- 

3 

modulated  sweep  utilized  in  the  P teclinique.  The  latter 
mode  employs  a sweep  of  pulses  which  follows  a sawtooth 
enveloi>e.  Between  pulses  tlie  sweep  returns  to  a constant 
baseline  voltage  Tlie  sequence  of  four  pulses  generates 

four  distinct  I-\  data  points  for  ttie  probe's  current- 
voltage  characteristic.  Ttie  fifth  pulse  is  blanked  out  so 
that  the  current  1 ^ collected  by  the  probe  during  the  inter- 
pulse period  can  be  monitored  and  used  as  a measure  of 
variations  in  the  probe-plasma  system.  These  variations 

might  include  clianges  in  the  probe-s  surface  condition, 

plasma  f luc:  t ua  t ions  , oi’  rapid  variations  in  plasma  density 

as  seen  fiy  a I'ocket  or  satellite  passing  through  an  in- 

homogenous  ionosi)heric  structure.  Tlie  values  for  T and 

on 

Tjj  (inset  lig.  2)  were  chosen  to  stabilize  probe  surface 
conditions  (Holmes  and  Szuszczewicz , 1975:  Szuszczewicz 
and  Holmes,  197.5)  and  establish  Ijj  as  a direct  measure  of 
plasma  density  variations.  Additionally,  the  ratio  of 
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reierence  electrode  area  (payload  plus  rocket  body)  to  that 
ot  the  sum  of  the  probe  and  its  guard  was  sufficiently 

4 

large  (~2(10  ) ) to  guarantee  that  area  effects 
( Szuszczewicz , 1972)  would  not  perturb  the  measurement 
procedure.  The  arrangement  of  electrical  parameters  in 
the  5-pulse  format  allowed  the  sampling  of  baseline 
current  Ig  every  9,75  msec  while  a complete  I-V  character- 
istic was  generated  in  380  msec. 

Analysis  procedure.  For  a probe  operating  in  the 
collisionless  regime  with  applied  electrical  potentials 
less  than  that  of  the  ambient  plasma,  the  current-density 
response  can  be  written  as 


IkT 

= N e / ® 

exp  (x„)  - J 

(-1 

J 2ffm 

^ 'J 

M. 

V e 

V 

1 

= .ie  ^ ii 

(la) 

(lb) 


where  Xp  normalized  difference  between  the  probe 

and  plasma  potentials,  and  respectively;  that  is 


X = e(V  -V  )/kT  , (2) 

'p  p o e 

P is  the  probe  radius  divided  by  the  electron  Debye 

length  Xg,  T is  the  ratio  of  ion-to-electron  temperature 

T./T  , m is  the  mass  of  an  electron,  M.  is  the  charged- 
1 e e 1 

normalized  ion  mass  M.  = m./Z®,  and  I.  is  the  dimension- 

1 i 1 


less  ion  current  defined  by 

.ii  - N^^e  x/kTg/27rm^  I^, 


(3) 
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In  Eqs.  (la)  and  (3)  tlie  quantities  as  yet  undefined 
are  the  magnitude  of  ctiarge  of  an  electron  e,  the  Boltzmann 
constant  k,  the  experimentally  observed  ion-current  density 
collected  by  the  probe  ,j  ^ , and  the  probe  area  a which 
collects  "effective  ion  currents"  as  a result  of  bulk 
plasma  flow  or  photoelectron  emission.  If  plasma  flow  or 
photoemission  effects  are  insignificant,  I^  can  be  taken 
from  the  woi'k  of  Laframboise  (1966),  otherwise  a dimension- 
less form  of  the  results  of  Wharton  and  Hoegy  (1971)  can  be 
applied . 

If  plasma  parameters  are  constant,  standard  analysis 

procedures  (Chen,  1965:  Boyd,  1968)  can  be  applied  to 

Eq . (la)  in  order  to  determine  tlie  electron  temperature 

and  density,  I'lie  procedure  involves  subtraction  of  the  j ^ 

component  I rom  .]  to  empirically  detei’mine  ,j  ^(V  ),  the 

e p 

first  term  on  the  right  in  Eq . (1).  Electron  temperature 

is  determined  from  the  slope  of  the  ln(j  ) vs  V nlot 

e p 

according  to  the  re  la  t ionsli  i]) 

r ->-1 


T 

e 


e 

k 


6 (In  ,j^) 
■ ^ I’ 


(-1) 


for  probe  potentials  greater  than  tliat  of  tlie  plasma 
(Xp>0)  electron  current  will  dominate,  and  the  current 
density  to  cylindrical  electrodes  can  be  represented  by 
(Chen,  1965) 

,j^(\  >0)“  N c 
e p V V 2iTm^  n j 
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1 


where 


k(\, 


h - 'll-  c 


xp  (\) 


1 - erf(xp 


and 


1 

\ f 

erf(X  2)  = 2 ^ gxp  (-t^ ) dt 

p Y ''  ' o 


(6) 


(7) 


For  Xp  ^ 4 and  Xp  ^ 10>  EQ • (5)  may  be  written  as 


Je  (Xp^4)  H 


'2N  e 
e 


yJ(kT^/2rrm^)  (1  + x ) 


(8) 


and 


Je(Xp^lO)  = yj2/ir^m^  ^ 

respectively.  An  identical  procedure  can  be  carried  out 

for  spherical  probes  with  appropriate  relationships 
substituted  for  Eqs,  (5),  (8)  and  (9).  Experimentally, 

the  electron  density  is  determined  by  the  slope  of  the 

appropriate  portion  of  the  I-V  characteristic  that  applies 

to  either  Eq . (8)  or  (9).  The  latter  has  the  advantage 

of  T^-independence . 

If  the  electron  density  varies,  T and  N can  be 

e e 

3 

determined  in  the  P technique  by  unfolding  the  density 
fluctuations  from  Eqs.  (1),  (4),  (8)  and  (9).  This  is  done 

with  inlerpulse  measurements  of  1^^  carried  out  at  a fixed 
potential  while  the  net  I-V  characteristic  is  collected 
along  the  sawtooth  envelope  shown  in  Fig.  2B.  Selecting 
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one 


so  that  I,.  IS  insensitive  to  T (e.g,  at  v ^ 1)  , 
[t  H e ’ 


Density-normalized  probe  characteristics  can  then  be  re- 

1*0  f 

constructed  from  tlie  raw  data  by  substituting  for 

everywhei'e  in  Kqs . (1)  througli  (9),  with 


N N ( 1-  , t ) B 

e e = ^ — 

I[j(r,  t)  . (11) 

In  practice,  this  procedure  involves  the  multiplication  of 

raw  cui-rent  values  on  the  characteristic  by  the  ratio 
1’  o 1 

^B  result  is  a curren  t-vol  tage  character- 

ref 

istic  i-efei’onced  to  a constant  plasma  density  . The 

selection  of  is  arl)itrary  but  generally  takes  on  one 

of  the  l|j-values  collected  during  the  generation  of  the 

I-V  characteristic  being  analyzed. 

ft  should  t)e  ix)inted  out  that  the  plasma  current 

collected  by  a probe  has  a number  of  subtle  responses  to 

density  vai’iations  wtiicti  were  not  discussed  abt)ve.  One  of 

t tu.'se  subtleties  is  the  eflect  of  density  on  the  value  of 

P and  its  associated  influence  on  I..  This  effect  is  noted 

1 

here  onlv  fOr  puiposes  of  completeness  because  tlie  pi’obe's 
(?-va  r ia  t i ons  ttiroughout  tlie  entire  flight  were  in  the 
or  bi  ta  1-mo  t ion- 1 imi  t ed  regime  (('hen,  1965;  Laframboise, 


196B)  and  tlierelore  did  not  alter  the  current  response. 
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This  was  not  necessarily  the  case  for  the  reference  elec- 
trode (the  payload  and  rocket)  which  operated  in  the  regime 
10  ~ /3  - 100,  Reference  electrode  variations  in  ^ would 
change  its  collected  ion  current  and  ultimately  its 
floating  potential.  It  has  been  shown  (Szuszczewicz , 1972) 
that  such  a ^-variation  might  result  in  a change  in  the 
vehicle's  floating  potential  equivalent  to  ~ 0,4kT^/e  volts. 
Variations  of  this  nature  were  not  present  in  the  sporadic- 
E application  reported  here, 

III-  RESULTS 

The  payload  was  launched  from  White  Sands  Missile 

Range  (32°30'N,  106° 20 'W),  New  Mexico,  on  29  June  1974, 

at  2031  MST  (0331  GMT,  30  June  1974)  and  reached  an  apogee 

of  229,5  km.  The  probe  was  deployed  upleg  at  102  km  and 
last  observed  ionospheric  plasma  on  reentry  at  90  km. 

Figure  1 shows  the  up-  and  downleg  profiles  of  relative 
electron  density  derived  from  the  probe's  saturation 
electron  current  collected  at  the  maximum  positive  sweep 
voltage  V^,  These  profiles  locate  the  up-  and  downleg  E^- 
layers  at  106,7  and  106,1  km,  respectively.  Within  the 
layers  the  probe  operated  in  the  pulse-modulated  mode 
(Fig.  2B) , providing  each  complete  I-V  characteristic  with 
running  measurements  of  relative  plasma  density  every  14 
meters , 
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i- i^;iire  3 is  an  expanded  view  ol  P currents  and  volt- 
ages as  I lie  payload  traversed  the  E^-layer  on  the  uple^t 

portion  ol  the  tli^ht.  Plotted  as  a lunction  of  altitude 

i c 

are  the  net  ion-  I and  electron  I currents  collected  by 

P P 

the  probe  during  voltage  pulses  which  followed  the  sawtooth 

envelope  (Fig.  2b).  The  baseline  current  values  I„  are 

also  shown  and  repj-esent  samples  of  relative  plasma  density. 

Since  the  probe's  iloating  potential  was  approximately  1 

volt  positive  with  respect  to  that  of  the  payload  (the 

reterence  electrode)  I samples  at  volts  were  collected 

in  the  ion-saturation  portion  of  the  1-V  characteristic  at 

a potential  approximately  equivalent  to  -25kT^'e  volts 

relative  to  the  ambient  plasma  potential. 

Because  plasma  density  is  a rapidly  changing  function 

ot  altitude  within  the  layers,  the  currents  and  were 

dei)endent  upon  the  local  density  as  well  as  the  voltage 
applied  to  the  probe  With  Ij^  providing  a measurement 

ol  variations  in  N , the  1 (V)  characteristic  was  unfolded 

e e 

1 rom  the  1 (V.N  ) data  bv  the  metliod  outlined  in  the  pre- 
e e ■ 

vious  section.  Then  the  ambient  electrui  temperature  T^ 

was  determined  in  the  I'e  t arc!  i ng-1  ield  region  of  the  probe's 

rc  1 

cha  rac  t er  i.s  t ic-  (see  Fq  . (1)  ).  For  each  1-V  curve,  1^^ 

wa.s  taken  to  be  that  value  ol  1 closest  (in  time)  to  the 
point  where  I*/  The  value  determined  for  T was 

as<i'ibed  to  that  |•eferenced  p(»sition,  that  is,  the  time 
and  altitude  ol  the  Iloating  potential  (defined  by  1^=1^). 
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Results  on  the  upleg  trajectory  yielded  (T^,Z)  = 
(340°K,  106.4km),  (345°K,  106,7km),  380°k,  107,5km)  and 
(515° K,  107.9km)  with  the  last  measurement  made  outside  the 
highest  density  region  of  the  layer  shown  in  Fig.  3. 

Identical  procedures  were  followed  on  the  downleg 
passage  through  an  layer  with  the  expanded-sca le 
currents  and  voltages  shown  in  Fig,  4.  Only  two  downleg 
measurements  were  made,  yielding  (T^,Z)  = (345° K, 

105.75km)  and  (500° K,  106,52km).  The  values  of  tempera- 
ture and  altitude  shown  within  parentheses  in  Fig,  4 are 
those  measured  in  the  upleg  E^-layer  and  superimposed  on 
the  downleg  profile  by  normalizing  to  identical  ratios  of 
( Z)/N^ (peak ) . Because  of  limitations  imposed  by 
electrometer  sensitivity,  reliable  measurements  of  elec- 
tron temperature  were  not  made  in  the  trough  region 

(109km  Z ^ 200  km)  above  the  Eg-layers.  Above  200  km, 
values  of  electron  temperature  were  between  850°  and  1000° K, 

Within  the  layers  absolute  values  for  N^(Z)  were 
determinable  from  the  slope  of  the  density-normalized 
saturation  portion  of  the  I-V  characteristics  appropriate 
to  Eq . (9),  This  analytical  technique  yielded  electron 
densities  of  3(10^)  and  7(10^)  cm  ® at  the  peaks  of  the  up- 
and  downleg  layers,  respectively.  All  other  values  of 
within  the  layers  can  be  scaled  from  the  Ig-profiles  in 
Figures  3 and  4. 
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The  diflerence  in  layer  thicknesses  (1.16km  upleg, 
0.73km  downleg)  is  in  agreement  with  the  generally  sys- 


I 


i 


tematic  downwai'd  motion  of  K^-layers  accompanied  by  the 
profile  growing  thinner  and  more  intense  (Layzer,  1972). 

If  in  fact  the  up-  and  downleg  observations  were  of  the 
same  layer  (the  horizontal  separation  between  up-  and  down- 
leg  Eg-layer  measurements  was  60  km) , the  average  downward 
drift  velocity  was  equivalent  to  1.8  m/sec. 

The  major  feature  of  our  observations  is  the 
measurement  of  rather  steep  electron  temperature  gradients 
across  the  layers  shown  in  Figures  3 and  4.  The  results 
show  the  bottom  side  of  the  layer  as  isothermal  at 
T^  as  345°  K while  the  topside  is  accompanied  by  a positive 
gradient  in  T^ , resulting  in  an  electron  temperature  of 
515° K at  the  top  edge  of  the  layer.  This  is  believed  to 
be  the  very  first  detailed  observation  of  electron  tempera- 
ture gradients  within  temperate  latitude  sporadic-E 

layers.  Previous  attempts  at  direc’t  T determinations 

e 

(Smith.  1966:  Aul)rev,  et  al,  1966)  whithin  E -lavers 

s 

were  not  as  detailed  as  the  results  waicti  we  present  here. 


IV.  DISCUSSION  OF  RESULTS 

a.  Coupling  to  Gravity  Waves 

In  this  section  we  consider  the  implications  of  our 
results  from  two  points  of  view.  The  first  assumes  total 
thermodynamic  equilibrium  in  the  charged  and  neutral 
particle  populations.  There  is  support  for  this  position 
(T  =T  =T  ) in  incoherent  radar  scatter  investigations  of 
the  daytime  ionosphere  in  the  undisturbed  region  (no  E^) 
below  130  km  (Carru  et  al.  1967a, b-  Evans,  1967-  Wand  and 
Perkins,  1968,  1970:  Wand,  1970:  Salah  et  al,  1975). 

Assuming  that  complete  thermodynamic  equilibrium  is  main- 
tained, electron  temperature  variations  are  traceable  to 
identical  variations  in  the  neutral  species. 

Variations  in  neutral  temperature  can  arise  from  the 

semiduirnal  tide  or  shorter  period  (<  4 hrs)  gravity  waves 

which  presumably  are  randomly  phased  fi’om  day  to  day 

(Harper  et  al,  1976-  Wand,  1976-  Salah  et  al,  1975  Hooke 

1969-  Hines,  1965).  The  semiduirnal  tide  however  is  an 

unlikely  candidate  to  explain  our  topside  E^  observations 

since  tidal  gradients  below  115  km  are  typically  in  the 

range  10°  km'*dT^^  clz^20°  km,  a value  hardly  comparable  to 

the  measurement,  dT  dz>30()°  Km. 

e 

If  one  is  to  ascribe  the  observed  T variation  to  the 

e 

influence  of  a changing  neutral  particle  energy  distribution, 
the  more  plausible  source  is  tlie  shorter  period  gravity 
wave  which  is  identified  with  the  wind  shears  giving  rise 
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to  sporadic-K  layers  (Hines  and  Co  1 leajiues , 1974).  Even 

here,  however,  there  are  protjlems  in  firmly  establishing 

the  connection,  particularly  t)ecause  the  data  base  on  the 

shorter  period  gi’avity  waves  i.s  virtually  non-existent  at 

night  and  undei’  conditions  of  sporadic-E. 

Typically  the  mid-latitude  incoherent  scatter  radar 

observations  of  1.  (assumed  equal  to  T ) have  been 

1 n 

confined  to  the  daytime  ionosphere  because  of  lack  of 
sensitivity  to  the  reduced  electron  densities.  Addition- 
ally, the  radar  measurements  are  influenced  by  T.,  T T.. 

. (ion-neutral  collision  frequency),  and  the  relative 
in 

concentration  and  masses  m^  of  ionic  species.  In  the 

referenced  radar  observa  t itins  a nominal  <m . > = 30.5  a.m.u. 

1 

was  assumed.  This  assumption  leads  to  a 2'<  error  in  the 

Tj^  determination  if  the  molecuilar  ion  NO*  (30  amu)  is  the 

only  heavy  charged  species  present  (Salah  et  al.  1975). 

This  sensitivity  of  the  radai'  measurements  to  the  m^- 

distribution  is  particularly  critical  under  E conditions 

s 

where  the  dominant  ions  can  be  Fe*(5(i  amu)  and  Mg^(24  amu). 

Another  difficulty  in  establishing  experimental 
support  for  a direct  cause-effect  relationship  between 
gravity  waves  and  our  1'^^  observations  i.s  that  our  data 
would  imply  a vertical  wavelength  ~ 1 km.  This  is  outside 
the  detection  capabilities  of  radars  which  typically 
integrate  over  a 2. 5- 3.0  km  range  and  require  0.5-1. 5 hrs. 
for  a complete  height  profile  of  temperature.  As  no  direct 
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measurements  of  fjravity  wave  parameters  were  made  during 
our  experiment,  we  can  only  speculate  on  the  characteristics 
of  such  a wave  at  the  time  of  our  observations. 

b.  Influences  of  Electron  and  Ion  Distributions 
In  contrast  to  the  proposition  of  total  thermodynamic 
equilibrium  and  heating  by  gravity  waves,  we  discuss  here 
the  effects  of  ion  composition  and  plasma  density  on 
electron  temperature  within  an  layer.  We  present  the 
predictions  of  Gleeson  and  Axfoid  (1967)  and  Hooke  (1969), 
extend  their  calculations  to  a multiple  ion-component 
plasma,  and  compare  the  results  with  our  experimental 
observa  t ions . 

Following  the  work  of  Gleeson  and  Axford  (1967)  it 
can  be  shown  that  an  E^-layer  with  two  ionic  species  ni 
and  n«  (and  associated  recombination  coefficients  tti  and 
dp),  has  an  altitude-dependent  temperature  differential 
between  electrons  and  neutrals  given  by 


-kT^  + 

o 

+ - Oa  ) 

di 

N 5 

<r>  - 2 

o 

kT 

e 

-1 

9^  Lg^(aiNo)"' 

is  defined  as  the  energetic  photoelectron  "injection 

energy"  (taken  by  Gleeson  and  Axford,  and  Hooke  to  be 

8000° K),  k is  the  Boltzmann  constant,  T is  the  neutral 

n 

particle  temperature,  is  the  ionospheric  plasma 
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density  in  the  absence  of  windshear.  and 


L*  = 1.55(10  "^)N(Ns)  4 [5(10“^®)T  - (10  ^^)]N(0b) 

^ (13) 

is  the  e lec t ron- 1 o-neu t ra  1 loss  rate  in  ev'  per  second  per  °K 
with  the  densities  of  Ns  and  given  in  number  per  cm^ . 

In  establishing  Eq . (12)  it  was  assumed  that  »ag  and 

(ni  )c^^(ns)~  , that  is  n;  »ns  in  the  absence  of  windshear. 
These  assumptions  are  supported  by  experimental  observations 
(Johnson,  1968;  Narcisi.  1968)  which  show  that  N, (=(ni), ) 
is  primarily  a molecular  ion  combination  of  mean  mass  near 
30  amu  (e.g.,  NO*  and  Oa ^ ) with  a recombination  coefficient 


4.2(10  ) (300  T^) 


Ui  ^ 7.5(10’')  (300  T^)cm^sec-^ 


(Oppenheimer  and  Brace,  1976  Oppenheimer,  et  al,  1976; 
Torr,  et  al,  1976).  i'he  second  ionic  species  n^  is  atomic, 
observed  witfiin  E^  layer.s  as  a metallic  (primarily  Fe*  and 
Mg^)  with  a very  small  radiative  recombination  rate.  ~ 
3(10~^^)  cm^sec”^  (Bioncli,  1972). 


If  only  molecular  ions  are  present  Eq.  (12)  becomes 

® 5 

2 


2kT' 
2 e 


‘n?’ 


kT 


(14) 


® " L*  (aNo  )~^  (N  /No  ) 

en  e 

Figure  5 presents  the  solutions  to  Eq . (14)  for 

altitudes  100,  115,  and  130  km,  “ 200°  K and  400°  K, 
aNo  - 3(10"^)  sec"^  and  (3  2)kT^=leV. 
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Table  1.  ALTilUI'i;  I)  I .''I  i:  I iU’ 1 1 ONS  ()l  MObKC'UI«M{  NITROGEN 


i 

1 A 1 1 i 1 udo  [ Km 

N'a  Dims  i t y [ cm  ] 

1 

I I 

I 1 I 

i 100 

2(  lo'  • ) 

2. 4 ( lo’  ‘ )* 

^2.5(10^^ )* 

1 lb 

1 . 4 ( lo’  ’ ) 

<;.2(10^^  )* 

^3.4(10^® ) 

130 

9 . 5 ( lo’  ° ) 

1.6(10^'  ) 

9.4(10^°) 

The  calculations  assume  the  molecular  nitrogen  densities 
shown  in  column  I ol  Table  1:  these  densities  result  from 
a fit  to  the  original  solutions  of  Gleeson  and  Axford  (1967) 
at  the  same  altitudes.  For  comparison,  the  Table  shows 
two  other  Na  distributions.  The  entries  in  column  II  have 
been  taken  from  Dalgarno,  et  al  (1967):  they  represent  a 
modified  Jacchia  1964  model  for  quiet,  undisturbed,  mid- 
latitude conditions  during  July  1963  at  1900  EST.  (When 
measured  in  terms  of  diurnal,  seasonal  and  solar-cycle 
influences,  these  conditions  identify  reasonably  well  with 
the  period  of  our  launch.)  In  column  III  are  the  rocket- 
borne  observations  of  Hedin,  et  al  (1964)  at  White  Sands 
Missile  Range  on  6 June  1963  at  0730  MST.  (The  asterisk 
in  columns  II  and  III  identify  extrapolated  values.) 

At  the  altitudes  in  question  we  have  neglected  the 
electron  energy  losses  to  molecular  oxygen.  With 
N(06  )<0.2N(Ns  ) and  T^^500°K,  Eq . (13)  shows  this  to  be  a 

reasonable  assumption  lor  simplifying  the  calculations  and 
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subsequent  discussions. 

The  minima  in  the  family  of  curves  shown  in  Fig,  5 
and  derived  from  Eq . (14)  occur  at  a critical  density 

ratio  given  approximately  by 


The  prediction  for  the  existence  of  a temperature  minimum 

within  an  E -layer  of  molecular  ions  means  that  as  one 
s 

moves  into  an  E^-layer  from  the  depleted  regions  outside 

(T  -T  ) should  decrease  as  N /N  increases.  If  the  layer 
e n e o 

were  sufficiently  intense  so  that  (N  /'N  ) , > (N  /N  ) , 

e o peak  e o c 

the  difference  (T  -T  ) would  achieve  a minimum  at  the 

e n 

critical  densitv  I’atio  and  then  begin  to  increase  as  N /N 

e o 

increased  toward  the  peak  of  the  layer.  Care  must  be 

exercised  in  evaluating  the  significance  of  this  prediction, 

particularly  since  high  recombination  losses  may  prevent 

the  development  of  a purely  molecular  ton  E^-layer  with 

(N  /N  ) , > 7,  A more  realistic  approach  includes 

e o peak 

metallic  ions  within  the  layer. 

To  demonstrate  the  effects  of  metallic  ions  on 
Eg-layer  temperatures  we  assume  an  ion  distribution  that 


conforms  to 


111  ■=  N,,  * 6(N  -N  ) 

^ ^ e o 


= 111  * ns 
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where  appropriate  selection  of  6 will  allow  either  an 
enhancement  (6>o)  or  depletion  (6<o)  of  the  molecular 
species  with  the  layer.  With  these  assumptions  Eq.  (12) 
becomes  _ 


T -T  = 
e n 


1 N 

1 N " 

1 

2 kT'  f 

1 . -1) 

1 _e_  5 

1 kT 

2 e 

1 

1 'TV 

1 N 2 

1 o 

e 

L*  (aiNo)~^  (N  /N  ) 
en  e o 


where  X = 6(l-cr^)  * arr 


Solutions  to  Eq.  (18)  are  shown  for  four  cases  in  Fig.  6 
where  the  running  parameter  X demonstrates  the  effects  of 
ion  distribution  (X<^6  for  realistic  values  of  tts/ai 
(«10“®  ) ) . 

The  case  X=1  represents  a purely  molecular  ion  layer 

where  a temperature  minimum  occurs  at  modest  values  of 

N /N  , As  X decreases,  the  contribution  of  metallic  ions 
e o 

to  the  total  plasma  density  increases,  and  the  location  of 
the  temperature  minimum  shifts  to  higher  values  of 
according  to  the  approximate  relationship 


Equation  (20)  applies  for  all  X>o,  For  X^o  there  is  no 
minimum. 

The  case  X=o  represents  the  transition  between 
molecular  enhancement  (X>o)  and  molecular  depletion 
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(X<o)  within  the  layer.  The  latter  case  is  perhaps  more 
representative  of  an  actual  condition,  being  in  agree- 
ment for  example  with  Young,  et  al  (1967)  who  found 
molecular  ion  concentrations  greatly  reduced  within  the 
layer.  In  contrast  to  the  case  of  the  pure  molecular- 
ion  layer,  this  suggests  that  predominately  metallic  layers 
with  molecular  species  depleted  by  enhanced  electron 
recombination  should  have  the  temperature  minimum  identified 
with  the  layer  peak. 

To  compare  the  predictions  of  theory  with  our  results, 

the  upleg-E^  profile  has  been  smoothed  and  presented  at 

the  bottom  of  Fig,  7 in  tlie  N /N  vs  Z format.  N has 

GO  o 

been  defined  by  the  I level  in  Fig.  3 in  the  region 

107 . 5szs 107 , 7 , This  delinition  and  the  analysis  procedure 

detailed  in  Section  2 yields  Nr,  =2 , 2 ( lo'^  ) cm~‘^  . 

The  difference  between  measui'cd  values  of  T and  an 

e 

assumed  neutral  temperature  of  300° K is  plotted  at  the  top 

of  Fig,  7,  with  open  circles  representing  upleg  data  and 

solid  circles  representing  downleg  values  normalized 

appropriately  to  The  error  bars  indicate  best 

estimates  of  Inaccuracy  in  the  T values.  The  estimate 

e 

stems  from  careful  scrutiny  of  Maxwellian  fits  to  retarded 
electron  currents  in  the  region  (\  -l)<xs;(x^,+  4).  The 
subscript  "f"  identifies  the  floating  potential.  The 
lowest  measured  value  of  T^  was  340° K and  the  largest 
estimated  error  was  ± 15° K,  With  the  procedure 
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stabilizing  probe  surlac:e  conditions  no  error  was 
identiiied  with  possible  measurement  distortion  by  con- 
tamination el  feels.  Tlie  selection  T^^=300°K  is  arbitrary, 
but  over  the  altitude  range  106  to  108  Km  it  is  consistent 
with  tiie  lower  t liermospher ic  temperature  results  reported 
by  Salali  and  Wand  (197'!). 

The  top  ol  Fig.  7 includes  the  predictions  of  Eq . (18) 

tor  the  parametric  cases  detailed  in  Table  2.  Curve  "a" 

IS  the  result  ol  applying  the  parameters  and  loss  terms 
used  by  Gleeson  and  Axford  (1967).  The  G-A  entry  in 
Fable  2 identifies  the  extrapolated  density  distribution 
which  follows  from  Table  1,  column  I. 

Table  2.  PARAMETERS  SELECTED  FOR  TWO-SPECIES  PREDICTIONS 

WITH  X-  0 AND  T =300° K 
n 


f “ 

I’a  I'amo  t er 

CURVE 

a 

b 

“^0  t J’GC  3 

3(10“^  ) 

1 .2(10“^ ) (300/T  ) 
e 

' [N,3 

G-A 

0.2  G-A 

' L* 

en 

Eq. (13) 

Eq.  (21) 

2kT'[eV] 
2 e 

. 

1 

2.0 
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Curve  "b"  represents  adjustments  in  the  parameters  in  "a". 
As  detailed  in  I’aljle  2.  tlie  adjustments  include: 

i)  The  use  of  a more  recent  analytic  representation 
for  electron  neutral  losses  (Banks  and  Kocharts,  1973; 
Dalgarno,  1969)  through  Ns  rotational  excitation.  This 
new  loss  term  is  given  by 

L*^  = 2.0(10^M  N(Ns)/.^^''  (21) 

with  an  uncertainity  of  a factor  of  2. 

il;  The  increase  of  the  electron  injection  energy 

skT^^  from  leV  to  2eV.  This  allows  the  upper  bound  on  T^ 

2 e e 

to  extend  beyond  the  l.SeV’  cut-off  identified  with  vibra- 
tional excitation  of  Na  in  the  energy  cascading  process  of 
photoelectrons  (Dalgarno,  et  al,  1963;  Dalgarno,  et  al, 
1967;  Timothy,  et  al,  1972).  When  one  includes  energy 
transfei’  to  thermal  electrons  from  vibra  t iona  1 ly  excited 
Na , the  greatest  upper  bound  on  T^  extends  as  high  as 
5eV  (Dalgarno,  et  al,  1963). 

iii)  A reduction  in  aj from  3(10  ^)  sec  ^ to 
1 , 2 ( lo”*"  )sec  ^ in  order  to  reflect  a value  more  in  keeping 

—7  3 —1 

with  experimental  c:onditions  (i.e.,  Qti  . 5 ( 10  ) cm  sec 

No  =2 . 2 ( lO”^  ) cm  ).  (!ase  ”1)"  also  reflects  temperature 
dependence  in  the  leaction  la  t e coeflicient  Ofi  . in 
accordance  with  Biondi  [1972],  Torr,  et  al  [1976],  and 
Oppenheimer  [1976]. 
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iv)  An  atljus  I men  t in  molecular  nitrot^en  densities 
whicli  l)riii!:.s  model  predictions  within  t lie  rant;e  of  our 
measu  remen  t s . 

In  evaluating;  the  theoretical  curves,  we  note  that 
ea<’h  could  liave  been  t;enerated  from  an  unlimited  number 
of  parametric  combinations.  Tliis  lack  of  uniqueness  in 
the  model  predic-tions  imposes  some  constraints  on 
conclusions  drawn  from  comparison  with  experimental  results 
For  example,  the  choices  for  and  Na  in  Case  "b"  were 

selected  to  improve  the  fit  between  theory  and  experiment. 
These  same  parameters  could  be  .justifiably  altered  to 
produce  solutions  that  ai’e  in  lesser  at;reement  with  our 
I’esults.  (We  note  that  the  Na  density  of  Case  "b"  repre- 
sents an  excursion  which  j;oes  beyond  existing  information 
on  neutral  density  variations  at  105  km  (See  e.g.  Faucher 
and  Morrissey.  1971).)  The  quantitative  agreement  is 
therefore  qualified  because  of  an  incomplete  determination 
of  all  existing  conditions. 

We  note  howevei'.  that  the  most  significant  feature 
of  (he  data  is  the  steep  gradient  in  electron  temperature. 
The  measured  t«-mperaturo  gradient  is  real,  and  its  relation 
ship  to  a G 1 eeson-Ax f ord  model  is  the  substance  of  Figure 
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SIMMAKY 


The  first  ionosplier  ie  a pi>  1 ica  t i<Mi  of  a newly-developed 
pu Ised-p lasma-prol)e  is  described  in  an  electron  temperature 
study  of  mid-latitude  sporadic-t.  The  major  new  result 
was  the  definitive  observation  of  an  electron  temperature 
gradient  across  the  layers.  Layer-peak  and  bottom-side 
values  of  T^^  were  approximately  ecjual  at  340°K  but  the 
topside  of  eacli  la  yer  stiowed  tempei  a tui  e increasing  to 
515°  K as  the  electron  <iensity  deci’eased  witfi  altitude  to 
the  undistui’bed  amfiient  t.-\egion  value. 

The  electron  t emptoa  I ure  ol>ser  va  t ions  were  compared 

witfi  t fie  ttieoretical  preii  ic  t ions  of  gravity  wave  tieating 

as  well  as  with  precficteif  efte<ls  of  variations  in  botfi 

electron  density  aiui  ion  composition  witfiin  t fie  sporadic- 

K layer  itself.  Witfiin  t fie  framework  of  plausifile  atmo- 

spfieric  variations  iieitfier  me<-fianism  could  fie  unimiely 

identifieef  witfi  t fie  ofjst'ivali  >ns.  Final  and  conclusive 

analyses  requirx'  simultaneity  in  t fie  measurements  of  , 

T , T.  , M..  and  T witfi  spatial  re.solution  muc;fi  greater 
ell  n 

than  currently  availafile. 
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FitJ.  2 — Continuous  :ind  modulated  sweep  modes  of  probe  opt'ration.  The  continuous, 
linear  sweep  represents  the  conventional  approach  while  the  modulated  sweep  schimiati- 
cally  presents  the  procedure  of  the  pulsed  plasma  probe  technique.  The  pulses 
were  applied  at  a 512  sec~'  rate,  allowing  for  19.1  pulses  during  the  voltage  sweep  be- 
tween V_  and  V.J.  The  inset  is  an  expanded  scale  of  two  adjacent  pulses. 
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Kin.  4 Kxpaiulcd  view  of  probe  vollant*  ami  currents  yieliimn  the  indicated  tem- 
perature measurements  of  34 5“ K and  50()'’K  at  tlie  bottom-  and  topside  of  the 
downlen  K,-layer.  'I'he  values  T,.  shown  in  parentheses  were  determineii  in  the 
uplen  layi’r  and  superimposed  on  the  downlen  profile  liy  identification  with  iden- 
tical ratios  of  . Values  of  Ijj  represent  the  density  profile  within  the 

layer. 


Fin.  <>  l)<‘|)cii(li'iK  t‘  of  li'mpcratim'  on  ion  i-omposilion  within  ;in  Kj,-layer. 
For  O'  • '<  1,  ^ +1  identifies  a purely  molecular  ion  layer.  iVs  X hecomes 

le.ss  than  +1  the  eontrihution  of  metallies  Im'couk's  inen'asmnly  more  impor- 
tiint  in  the  total  ion  population.  A < 0 identifies  metallic  layers  accom- 
panied hy  molecular  (let)lelion. 


